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An analysis of the energy spectrum and the magnetization curve of two-dimensional organic anti- 
| ferromagnet F2PNNNO with a spin-one dimerized structure shows that a behavior of the compound 

in an external magnetic field can be explained within a lattice boson model with an extended Pauli's 
exclusion principle, i.e. no more than two bosons per a dimer. The unusual magnetization curve 
observed experimentally in the compound reflects a sequence of phase transitions intrinsic for a 
lattice boson system with strong on-site and inter-site repulsions due to a tuning of magnon density 
■ by the applied magnetic field. 
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I. INTRODUCTION 



A possibility to study the Bose-Einstein condensation (BEC) with low-dimensional magnetic materials predicted 
theoretically twenty years ago! gave rise to intense experimental studies in the field. The analogy between the 
spins and the bosons becomes evident for antiferromagnets where spins form dimers with a spin-singlet ground state.— 
Originally, main attention was focused on spin-1/2 systems where excitations inside each dimer (triplons) are regarded 
as bosons with hard-core repulsion, i.e. no more than one boson present on a single dimer. The analogy enables to 
treat the spin systems as that of interacting bosons whose ground state is determined by the balance between the 
kinetic energy and the repulsive interactions. 3 If the repulsion dominates the bosons will form a superlattice and a 
finite energy cost is needed to create an additional particle. This exhibits itself as a jump in chemical potential versus 
boson number, in the spin language, as a plateau in magnetization curve versus magnetic field at rational fraction of 
saturated magnetization. 

I \ The field induced condensation of magnons has been experimentally observed in coupled quantum (s = 1/2) dimer 
■ systems based on Cu 2+ ions such as TlCuCbj and BaCuSi2 0g^££ and the compound Ba2Cr 2 0g£ which are adequately 
Q> described by the BEC theory. 

Recently, magnetic weakly coupled dimer system BaaM^Os with 5 = 1 moments attracted a lot of attention.— ^ The 
field behavior of magnetization in the system of antiferromagnetically weakly coupled 5 = 1 dimers can be described 
as BEC of magnons by mapping the spin-1 system into a gas of semi-hard-core bosons*^ On an example of simple 
t-H ■ two-dimensional (2D) 5 = 1 isotropic Heisenberg model with a dimerized structure and frustrating interactions it was 
suggested an emergence of the spin supersolid state (a long-range mixing of superfluid and charge ordered phases) 
induced by a magnetic fields 

The organic compound F2PNNNO is a supplementary example of spin-one dimer based magnetic insulator. This 
. is 2D Heisenberg system with a singlet ground state, in which 5=1 dimers interact antiferromagnetically ; 12 i 13 
The lattice of the system is equivalent to the honeycomb one (Fig. fj}. The field magnetization process shows 
a two-step saturation behavior that is a rare example of observation of a plateau in a two-dimensional system. 
The intermediate plateau corresponds to the half value of saturation magnetization. The consistent calculation 
of susceptibility and magnetization for the finite-size cluster with imposed periodic conditions yields the following 
estimations of antiferromagnetic exchange couplings 2 Jo = 67.5 K, 2 Ji = 7.5 K, i.e. the system can be regarded as a 
real 2D dimerized spin-one system. 

Apparently, the quantum antiferromagnet F2PNNNO offers an opportunity to verify a relevance of semi- hard core 
boson model for description of the dimerized system. In the paper we perform a diagonalization of finite cluster 
of N — 18 sites, calculate the magnetization and demonstrate that these results can be easily understood within 
the semi-hard boson model with strong on-site and inter-site repulsions. The diagonalization procedure used by us 
accounts the spin rotational symmetry ] 14 i 15 The implementation of non-Abelian SU(2) spin symmetry is based on an 
elimination of quantum numbers via the Wigner-Eckart theorem. The advantage of the approach is that the cluster 
spin states are decomposed into different sectors of the total cluster spin. In addition, one can independently handle 
each of the target spin state. 
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The paper is organized as follows. The model and the diagonalization algorithm are given in Sec. II. The truncation 
procedure is discussed in Sec. III. In Sec. IV we report numerical cluster calculations of the spectrum and the magne- 
tization curve. The analogy with the lattice boson model is performed in Sec.V. Main results are recapitulated in the 
Conclusion part. 



II. THE MODEL 



The Hamiltonian of weakly interacting spin-one dimers on a 2D lattice depicted in Fig. [T] is given by 

Hs = Jo / ] SjlSj2 + J\ 2_j SiaSja, (1) 
i (ia,ja) 

where Jo is the coupling inside the i-th dimer, J\ is the strength of the exchange interaction between the dimers 
located on the bonds The indices a, a mark S = 1 spins that enter into the interacting dimers, namely, a 

= 1,2 provided a — 2,1, respectively. The both types of the interactions are antiferromagnetic Jo,i > 0, and the 
regime of weakly interacting dimers, |«7o| 3> | Ji|, is considered. The Heisenberg model has been previously suggested 
to explain some thermodynamical properties of F2PNNNO. 12 Numerical calculations based on the Hamiltonian ^ 
via exact diagonalization of small clusters and their comparison with experimental data prove its relevance for the 
ratio I Ji/Jq\ <C 1. The dimerization caused by the anisotropy of interactions on a lattice is somewhat analogous to 
a situation in two-leg spin-1 antiferromagnetic ladders in a strong antiferromagnetic rung-coupling regime, when the 
ladder ground state is well approximated by the tensor product of singlet rung-dimersJ^ 

To get the energy spectrum the finite-size clusters composed of N — 10 and N — 18 sites are selected. In a 
choice of the cluster care should be taken to ensure that the lattice point group symmetry is hold. Since intra- 
dimer interactions are the strongest, the cluster should consist of intact dimers. To mark sites inside the cluster the 
chessboard-like notations will be used, i.e. site positions along the x axis are marked by numbers whereas positions 
along the y axis are denoted by Latin letters. 

To find eigenfunctions of the cluster that inherit the total cluster spin as a quantum number we should develop a 
consecutive procedure of addition of spin moments. It is convenient to break the cluster in several parts. Following 
the strategy of a cluster building used in RefJ^, one should identify the central dimer (center) and its environment. 
The center is composed of the c3 and d3 sites whereas another sites are embodied into the environment. 

The Hamiltonian of the central dimer has the form H c — JaS C 3Sd3, whereas the interaction between the center and 
its environment is given by 

V ce = JlS c 3 (§ c 2 + S c ij + JiSd3 (Sd2 + Sdij ■ (2) 

The environment consists of four parts, namely of two dimers, left (I) and right (r) ones, with the Hamiltonians 

Hi = JoSdSdi, and H r = J^S^Sd^, (3) 
respectively, as well as two fork-like parts, i.e. the down and upper ones, with the corresponding Hamiltonians 

ffdown = Jo I Sb2S c 2 + S a 3Sb3 + Su^ci ) + JlSb3 I Sb2 + 5*64 ) , (4) 



H up — Jq ySd2S e 2 + S e sSf3 + SdiSeij + JlS e 3 yS e 2 + S e ij . (5) 

The interaction between the left/right dimers and the fork-like parts is presented as 

Knv — J\ 

The Hamiltonian of the entire cluster gathers all the above terms 

H = H c + V cc + {Hi + H r + Hdown + ^up + V cnv } . (7) 

There are three states of the dimer, which is the elementary block of the cluster, with the total spin Sd m = 
(singlet), Sdm = 1 (triplet), and Sdm = 2 (quintiplet) . The energies of the states are Eq = — 2J , E\ = — Jo, E 2 = Jo, 
respectively, and the eigenstates are obtained via the common rule of addition of moments 



1 11; SdmMim) = l^dm-Wdm) = / J 



1 1 Sdm 
01 02 Mdu 



IW |ltr a ) , (8) 



1 2 3 4 5 



FIG. 1: The 18-site cluster used in numerical calculations. The environment of the central dimer consists of two "fork"-like 
parts (up and down), and the left (I) and right (r) dimers. The intra-dimer Jo and inter-dimer Ji interactions are shown by 
solid and dotted lines, respectively. 

where [. . .] is the Clebsch-Gordan coefficient. To increase the cluster size the reduced matrix elements (RME) of the 
spin operators 5(1) and 5(2), that constitute the dimer, calculated within the basis (|5J) are needed 

(5 dm ||5(1)|| 5 dm ) = (-1) 1+5 - yJ(2S dm + l)(25 dm + 1) | 5 f \ 5 ^|(1||5||1), (9) 

(S dm ||5(2)|| S' dm ) = (-l) 1+s - v /(25 rfm + l)(25^ + l)| Sd ™ J S *» } (1 1|5|| 1) , (10) 

where {. . .} is the 6j-symbol of the rotation group, and the reduced matrix element (1 ||5|j 1) = \/6. 

The fork-like part includes three interacting dimers. It is convenient to build the basis of the fragment according to 
the scheme (2 + 4) + 3 of the moment addition, i.e. a combining of the "prong" dimer functions is followed by adding 
of the "handle" function. As a result, the basic functions with the total spin 5d OW n of the down fork-like part has the 
form 

|(5254) 524, 53; 5d own A/down) = 

M2M3M4M24 

where 52, 53 and 54 are the spins of the dimers composed of the b2 and c2 sites, etc. Within the basis the Hamiltonian 
(H|) is presented by the block diagonal 141 x 141 matrix. The blocks are marked by the total spin 5down = 0, 1, . . . , 6 
values. A diagonalization of the -ffdown matrix yields the spectrum -Ei down s down and the eigenfunctions 

|*down5down-^down) = ^ ' Oi ( S2S4 ) ^24', S3 I (^2^4) $24) 5 Sdown A^down) , 

52 S3S4S24 

where the index «down distinguishes basic functions with the same total 5d OW n spin. The results for the upper fork-like 
part can be obtained analogously provided the site c4 is substituted for d2, and c2 is changed by dA etc. The assembly 
of the cluster part is completed by calculations of the reduced matrix elements [see Eq. (|All) in Appendix]. 



52 54 524 

M 2 A'h M 24 



524 53 5d 

M 24 M 3 M d 

own 



|5 2 Af 2 ) |5 3 Af 3 ) |5 4 Af 4 



(11) 
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As the next step, we construct the spin functions of the non-interacting parts, i.e. of the left and the right dimers 

\S1M1) \S r M r } , (12) 



\SiS r ; Sir Mi r ) = 

MiM r 



Si S r Sir 
Ml M r M W 



where Si r = 0, 1 . . . , 4, and the upper and down fork-like parts 



l^up^upMownSdowni S u dM ua \) — ^ ^ 

M up M do 



^up ^down S u d 

M U p M([ own M u( { 



l^upSup-^up) l^down'S'down-^down) , (1^) 

where S n d = 0, 1 . . . , 12, and add them together to build the basis of the environment of the central dimer 



I (^up'S'up^down'S'down) S uo \ : {SlS r ) Sl r: S env M^ e 



E 

M ud M ir 



Sud Sir S cnv 

M u( i M^ M cnv 



|*upS'upidown»5'down ! S u dM u d) | ^ 5" r ! Sl r Mi r ) 



(14) 



The reduced matrix elements of spin operators needed to build the Hamiltonian of the environment are relegated to 
Appendix [see Eqs. (|A2IIA5j) ]. Note, that a number of the states (fT4|) is too much to avoid the truncation procedure 
(see Sec.III). 

Matrix elements of the environment Hamiltonian H onv = Hi + H r + i^down + H np + V eav are given as follows 

{(^up^up^down^down) *Sud} (SlS r ) Sl r , S env M^ env \H env \ (^up^up^down^down) ^ud' (S/'-V) *^/r' ^env^env) 

= (-EiupSup + E idamxSdom + E Sl + E Sr ) Si^Js^s^Ji^j^ds^^s'^d 



TX I 1 \S eBV +S' A +Sir- J "^ud Sir Sen 

•AOSenv,.?' ("I) 1 Qf C" 1 

I D lr °ud 1 



Sm„, 



ML 



(15) 



X {(SlS r ' t Sir 1 1 ■S'c 1 1 1 SiS r ', Si r ) ^JupS'upidown'S'down; Sud ||»S C 2|| ^up^up^down'-'down' ^ud) 

+ (SiSr, Sir \\Sdi\\ SiS r ; S ir ) (iupS' U pidown*S'dowii; S u d \\Sds\\ * up 5' up i down iS' down ; S ud ) 

+ (SlS r ', Sl r 1 1 'S'c5 1 1 SiS r ', Si r ) ^iupSup^downSdown; S u d 1 1 *S'c4 1 1 * U p'-'up*down'-'down i ^ud) 



+ (SlS r \ Sir \\Sd5W S[S r \ S[ r ) (iup'5'up^down'S'down; S u d ||'S'd4|| *up'^up*down'-'down i ^ud) } ■ 



The terms in {. . .} include product of the reduced matrix elements given by Eqs. (IA2IA3|) for spins that enter into the 
left /right dimers and by Eq. (|A4IA5j) for the constituents of the fork- like parts. 
After a finding of the environment eigenvalues Ei ciiv s mv and eigenfunctions 



I ^ e n vSc n v M c n v ) 



JenvSanv 

(iupSupidownSdowii)S„d 



(SiSr)Si r I (^upSup^down^down) ^ud; {SlS r ) Sl r ', S cnv M cnv ) , (16) 



one calculate within the basis the reduced matrix elements for the environment spins that directly interact with the 
central dimer see [Eq. (|A6[) ]. 

At the final step of the diagonalization procedure one build the basis of the entire cluster 



sScnv, Sc, SM) — ^ ^ 



M cnv M c 



Scnv S c S 

M env M c M 



l^env'S'onv-^cnv) \S C M C ) , 



and determine the matrix elements of the cluster Hamiltonian ([7]) 

(ienv^env, S c ] SM\H\i' env S' env , S' c ; S' M') = (E icnvSonv + E S J S imv y Ss env ,s' Ss c ,S'Js,S'Smm' 



5 

T ( i\S+s'+s c J S cnv S c S \ x x 

+ Jl( — l) onv < Q ; Q , , >0s.S>OM,M> 
I °c °cnv 1 J 



<S C ||S(1)||#) < 

^cnvScnv \\Sk\\ * C nv^cnv) ~t~ 11^(2)11 S c ) ^ ] (ienv^Jeaav 'cnv^ci 

fc=d2,d4 



fc=c2,c4 



(17) 



where the RMEs are previously derived [see Eqs.([9l fl0|) and Eq. (|A6|l ]. Numerical diagonalization of the matrix (fT7|) 
yields the target spectrum Eis and the eigenfunctions 



\iSM) =^lf mMm ,a e l^c„vM cnv , S c ; SM) 



(18) 



III. TRUNCATION PROCEDURE 



The classification of eigenstates of parts used to gather the total cluster according to irreducible representations of 
S'J7(2)-group enables to organize a truncation procedure inside the sectors of Hilbert space that arise at the consecutive 
steps of the algorithm. A possibility to carry out calculations within a reduced basis is feature of the algorithm that 
relates it with other renormalization group methods. 

We hold the following strategy of the truncation procedure to build target states that are obtained after combining 
two parts of the lattice. For a given spin-5* sector a certain amount of states having the lowest energy are kept. Thus 
every group of the \iS) states is presented in the reduced basis. We truncate the basis of two "fork"-like parts before 
to combine them into a larger lattice segment. This is not the unique way, for example one can truncate the basis of 
the environment after combining the "fork" -like parts, but the former is easier to perform. 

We tested several realizations of the truncation procedure either by simply controlling a number of vectors retained 
in the reduced basis or by monitoring a genealogy of the target spin-5 1 state through the triangle rule, i.e. only states 
that contribute into the target state are took into account. The last approach gives an opportunity to keep more 
vectors in the basis due to an omitting of redundant states. Moreover, highest-spin cluster states, i.e. those with 
S > 15 in our problem, are treated exactly. The size of truncated basis was chosen to be equal to either 64 or 121 for 
the scheme without an account of genealogy of the target state, and it varies from 12 till 352, being depending on the 
total spin S, for the "genealogical" scheme. 

An accuracy of the truncation procedure is controlled by monitoring an energy of the lowest state within an each 
spin sector. The variation of this observable computed through the both schemes does not usually exceed 1-2% (a 
maximal discrepancy of order 6% is reached only in the S-8 sector) that evidences a correctness of the constructed 
basis which exhibits almost no dependence on the used truncation procedure. The results that we present below are 
obtained within the "genealogical" scheme. 

Another feature of the algorithm is an addition of a central unit (one site or dimer) with its environment at the final 
step. The procedure does not depend on a structure of the environment and looks similar for any cluster. However, 
the information about quantum numbers of the environment states enables to simplify calculations substantially at 
the stage of the algorithm. Indeed, for a given spin-S* sector of the Hilbert space of the entire cluster one should pick 
out only those environment eigenfunctions whose spins S u obey the rule 

\S U S c \ < S <• S u + S c . 

Using of the truncation procedure results in the basises composed maximally from 4-5 thousand states. To control 
an accuracy of the procedure the results obtained for the 18-site system are compared with those for 10-site system. 
The smaller cluster enables to handle a whole basis without any truncation. The 10-site system is embedded into 
the bigger cluster and consists of the following parts: the central dimer c3, d3 and the neighbor dimers b2, c2, &4, c4, 
d2, e2 and dA, e4. Apparently, a construction of the environment requires two consecutive steps (i) an addition of 
the dimers &2,c2 and b4, c4 as well as <i2,e2 and dA, e4 ones according to Eq. (fl"2"|) then followed by a calculation of 
reduced matrix elements according to Eq.(A2 A3); (ii) a construction of the environment, states from the upper and 
down parts built previously and a calculation of RME of the environment spins that interact directly with the central 
dimer. The entire cluster Hamiltonian is obtained through (fTT)i . The biggest Hilbert space dimension (2025 x 2025) 
is reached in the S-2 sector. The numerical results for the supplementary cluster are listed in Table I for comparison. 
Note that one should compare energy values with the same magnetization per dimer (See Fig. [2]). 
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FIG. 2: Plot of the lowest energy per dimer e (to) vs m for the N = 10 and N — 18 clusters. The cusp is seen at m = 1. 



IV. ENERGY SPECTRUM AND MAGNETIZATION CURVE 



The results of the energy spectrum calculation for two N = 10 and N = 18 clusters are listed in Table I, where 
we give minimal energy E m i n within the each spin-S 1 sector along with the energy per dimer e = 2E m \ n /N . The 
magnetization per dimer is determined by m — 2S/N. The N = 10 and N — 18 dependencies e (m) are shown 
together in Fig. [2] The points for both clusters lay on one curve, i.e. finite-size effects may be ignored that is 
expected for the regime of a small dimer-dimer interaction J\ <C Jo- 

A remarkable feature of the curve is a cusp in the middle, i.e. at m — 1. The independent fitting of both parts 
jointed in the point by the quadratic form e(m) = E2m 2 +eim + eo yields £2 = 0.190 ±0.018, ei = 0.828 ±0.019, and 
e = -2.0073±0.0040for the lower part of the curve (0 < m < 1) together with e 2 = 0.200±0.058, £1 = 1.4578±0.018, 
and £ = -2.629 ± 0.014 for the upper part (1 < m < 2). 

On the base of data for N=18 case we build a dependence of jumps E m i n when the total spin S changes from till 
18, or the dimer magnetization varies from till 2 (Fig. [3]) One can see that the values of the jumps are approximately 
Jo for S < 9 and they are increased by a factor of 2 as S > 10. It means that the energy of the total system of 
weakly interacting dimers will change with an increasing magnetic field due to local excitations inside separate dimers. 
Indeed, for the single S = 1 dimer the spectrum consists of a singlet, a triplet, and a quintuplet. The energy difference 
between the singlet and triplet is Jo while the difference between the quintuplet and the triplet is 2 Jo (see discussion 
in the next Section). 

A standard way to consider a magnetization process at T — is to define E^ n (N) as the lowest energy of the 
Hamiltonian ([1} in the spin-S" subspace for a finite system of N elementary dimers. Applying a magnetic field 

B leads to a Zeeman splitting of the energy levels E^^(B) = E^ in — SB, and therefore level crossing occurs at 

values Bs — E min (B) — E^ a {B) when increasing the field. These level crossings correspond to jumps in the 
magnetization at zero temperature of the value 1/iV, until the fully polarized state with the magnetization per dimer 

m sa t = 2N/N = 2 is reached at the value of the magnetic field B sett — E^^(B) — E^-^ X \b). The calculation 
performed for N/2 — 9 dimers yields the magnetization points presented in Fig. 0]and demonstrates an appearance 
of the ground state plateau as well as the plateau at one-half of the saturation value. 

To guarantee a validity of the magnetization curve we use an approach developed by Sakai and Tahakashiii to recover 
the m(B) dependence in the thermodynamical limit. In this case the condition for the crossover fields transforms into 
B = e'(m), where e is the energy per dimer. The plateau boundaries are determined by the derivatives in the special 
points: (i) B\ — £'(+0) is related with the end of the ground state plateau; (ii) Bi = £'(1 — 0) and B3 = e'(l + 0) 
correspond to the beginning and the end of the intermediate plateau, respectively; (iii) B4 = e'{2 — 0) marks an 
emergence of the saturation magnetization. 

A treatment of the energy spectrum results in the linear dependences relevant for sectors between the plateaus 

e'(m) = 0.83 + 0.38 m, < m < 1, , . 

e'(m) = 1.46 + 0.40 m, 1 < m < 2, 
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FIG. 3: Plot of the changes A.E m j n versus the dimer magnetization m. A distinct jump is seen at m = 1. 



TABLE I: Numerical results of the lowest energy E m - m and the energy e per dimer in the spin-5* subspaces for N = 10 and 
N = 18 clusters. 



s 


E min {N = 10)/J 


i(N= 10) 


£ min (iV = 18)/Jb 


e(N= 18) 





-10.0334 


-2.0067 


-18.0336 


-2.0037 


1 


-9.1853 


-1.8371 


-17.1431 


-1.9048 


2 


-8.2123 


-1.6425 


-16.2529 


-1.8059 


3 


-7.1978 


-1.4396 


-15.2935 


-1.6993 


4 


-6.1430 


-1.2286 


-14.3205 


-1.5912 


5 


-4.9344 


-0.9869 


-13.3164 


-1.4796 


6 


-2.9787 


-0.5957 


-12.2745 


-1.3638 


7 


-0.9610 


-0.1922 


-11.1879 


-1.2431 


8 


1.0849 


0.2170 


-10.0260 


-1.1140 


9 


3.1588 


0.6318 


-8.8335 


-0.9815 


10 


5.4418 


1.0883 


-6.8807 


-0.7645 


11 






-4.8994 


-0.5444 


12 






-2.8795 


-0.3199 


13 






-0.8172 


-0.0908 


14 






1.2815 


0.1424 


15 






3.4533 


0.3837 


16 






5.6844 


0.6316 


17 






7.960 


0.8844 


18 






10.3254 


1.1473 



that produces immediately B± = 0.83 Jo, B2 = 1.21 Jo, B3 = 1.86 Jo, and B4 = 2.26 Jo- The values normalized to 
the saturation field B sat are listed in the Table II and exhibit a reasonable agreement with the experimental data for 
F2PNNNO system. A comparison of the finite cluster calculations with those of the thermodynamical limit (| 19|) is 
given in Fig. 2J It is seen that both methods produce the close results. 

Note that the method we used for numerical calculations is intrinsically two-dimensional one whereas the previous 
numerical study of the system^ dealt with the cluster embedded into a chain. The regions between the plateaus of 
the magnetization curve exhibit a behavior closer to linear one instead of the S-shape forms early obtained. 
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FIG. 4: Plot of m versus B obtained via B — e'(m). The dots mark values found through the diagonalization algorithm. 



TABLE II: Values of the magnetic field special points compared with the experimental data. 



Bi 1 B sat 


i = 1 


i = 2 


i = 3 


i = 4 


Theory 


0.37 


0.53 


0.82 


1 


Experiment^ 


0.33 


0.53 


0.89 


1 



SEMI-HARD CORE BOSON MODEL 



Let us introduce the boson picture on the base of data presented in Fig. [3J For J\ <C Jo the low energy subspace 
of spin Hamiltonian (fTJ) consists of the singlet, the S z — 1 component of the triplet, and the S z ~ 2 component 
of the quintuplet. It is convenient to identify the triplet state with a presence of a bosonic particle (triplon), the 
quintiplct state as a pair of bosons (quintuplon), and the singlet state as an absence of bosons. Then the boson model 
is formulated via the semi-hard core bosonic operators <?; and g\ with the extended Pauli's exclusion principle g\ ' = 0, 
i.e. more then two bosons per site are forbidden. Note that the principle may be realized via parafermion language 
but the description requires a transmutation of statistics that complicates calculations in a 2D case (see Appendix 



B). The algebra of the operators are [gi,gi 



9i,9i 



= 0, and 



— 5ij (1 — Fi), where Fi = (3/2) rii (rii — 1) is 



the deformation of the canonical boson algebra, rii = g\ 9i is the number operator 
The boson Hamiltonian in terms of these operators is written as 

U 
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H 



2 ^ 

(ij) 



9i9j 



9)9, 



(hi 



^2n, + — ^2n % ( ni + 1) 



{ij) 



1) K 



(20) 



where the hopping terms 

hi = t\ (n l j 



2) (riij - 3) , h 2 = 2t 2 {riij - 1) (3 - n tj ) , hi = t 3 (n tj - 1) {m 



depend on a number of particles rii 



rij on the bonds 



5) 



The map between the bosonic (T2TH) an d the spin Hamiltonian (fTJ) is reached through the representation^ (see Fig. 



5Z QZ 1 QZ 

i — °il "r D i: 



where (1,2) marks two spins on each dimer, and 



r c rS 



n=0.5 



2Jn 



S=l 



Jo 



s=o- 



CO-phase 



N \ 
\ 

\ 

\ 

\ 



.V. 

t r 

f Mi-phase 



n=1.0 



,S Z =2 

S=0 
■ S=l 



FIG. 5: The low energy subspace of the single dimer spectrum in the presence of a magnetic field. Boson superlattice patterns 
corresponding to the charge-ordered and Mott insulating phases are shown above. 



This establishes the relationship between the spin and the bosonic parameters U = Jq, V = Ji/2, /i = B — 4Ji, and 



U 



8\/2 
" 3^ 



a 1 Ji, where a - 



1, 2, 3). Thus, the bosonic model includes the strong on-site boson repulsion U as 



well as the noticeable repulsive intersite interaction V. The magnetic field B plays the role of the chemical potential 
fi. 

The boson Hamiltonian (f2T)|) constitutes low-energy effective model of the spin Hamiltonian |T]) that appears from 
restricting Hs to the subspace of the semi- hard core bosonic operators. The map is valid in the limit J\ <C Jo, or 
in the boson language ti/U, V/U <C 1, when the main physics is governed by a competition between the one-site 
repulsion and the chemical potential. 

The quantum phase diagram of the boson Hamiltonian (|20|) was built in Ref.— by using the stochastic series 
expansion quantum Monte Carlo method (see there Fig. 4). It has been found that a Bose condensate fraction appears 
in the regions of the chemical potential (magnetic field) between the platos of the <?-particles density (magnetization 
curve). In contrast, the charge density wave (Ising-like charge order (CO) phase) forms around the intermediate 
plato. There are regions, where supersolid phase, a mixing of the charge order and the Bose-supefluid (BS), emerges. 
According to the study the magnetization curve shown in Fig. 4 can be interpreted as a tuning of boson density 
by the applied magnetic field. At small chemical potential the empty states has the lowest energy, when all dimers 
in the singlet state (boson vacuum). For B > B\ a finite density of bosons (triplons) emerges in the ground state 
and contributes into a BS phase. The triplon excitations are mobile due to weak interdimer coupling. The density 
(magnetization) increases monotonically as a function of magnetic field until B2, where a transition to the CO-phase 
comes up. This corresponds to the boson concentration n — 0.5, when the triplons crystallize in a superstructure 
pattern (Fig. [5]). The fractional plateau requires strong boson interactions in comparison to the kinetic energy. At 
B > B3 the filling increases monotonically in the resulting BS phase (quintiplon condensation) till the ground state 
transforms into a Mott insulating (MI) phase with two bosons per dimer at B > B±. The boson concentration in the 
MI phase n = 1. The reasonings are easily reproduced if to analyze the boson Hamiltonian ([2H|) by neglecting the 
intersite terms. 



VI. CONCLUSIONS 



Quantum dimer antiferromagnetic systems is a nice testing area to study BEC of interacting particles. Along with 
the ultracold atomic gases in optical lattice o 18 ' 19 they offer an opportunity to observe transitions predicted by lattice 
boson models. In many problems the boson picture is more physically transparent than the original spin language. On 
the base of the analysis of the finite cluster energy spectrum for the two-dimensional spin-1 organic antiferromagnet 
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F2PNNNO with the dimerized structure we prove a relevance of the model of semi-hard core bosons with pronounced 
on-site and inter-site repulsions for the low-dimensional spin system. The unusual magnetization curve observed in 
F2PNNNO is nothing but a manifestation of fine-tune of density of the bosons by the applied magnetic field, when 
the low-density Bose-superfluid, charge ordering with one boson per a dimer, and high-density Bose-supefluid phases 
change subsequently each other with an increasing of the field. 
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APPENDIX A 

The reduced matrix elements for spins on the c2 and c4 sites computed in the basis of eigenfunctions of the 
Hamiltonian i?down are given by the 141 x 141 matrix 



idown'S'down || <5*c2(c4) || *down'^'down / — ^ ^ n .'c , c ^ '1 



''downSdown „'do™ S down 

j "(S 2 S 4 )S 24 ,S 3 a (S^i)S^,5^ 
i en •> ' 



X ((S2S4) S24, S3, Sdown || <S'c2(e4) || (S2S4) ^24; ^3) ^dovm) ■ (Al) 

The reduced matrix elements that enter into the expression are calculated according to the rules 

((S2S4) S24, S3; 5down \\Scs\\ (S2S4) S 2 4, S' 3 ; S' down ) = (— i) S2 + S4 + S3 + S24 + s, 24+ 5 'down [5 24; g' 2 ^ SdoWIl, S'downl ^ 

X \ S «> I s ' 2 A\ Sd °r I s '^\(n-,s 2 \\s(2)\\n-,s! 2 )Ss iSi Ss 3 s' 3 , 

[ ^2 A l ^2 J [ "J24 ^3 *J24 J 
((S2S4) S24, S3; S'down ||<Sc4|| (^2^4) ^24> ^3) "^dowil) = (— l) S2+ ' S4+ ' S3+2S,24 + S,down [S'24, 5 2 4, S'down, S down ] ^ 

x{ 5 Q 2 , 4 I S ?)\ S %r I ^)<H;S4||S(2)||li;5i>Ww 

^ D 4 t>2 04 J [ *24 *3 <J24 J 

where [S] = (25+1). 

The reduced matrix elements for spins on the sites cl (til) are given by the 19 x 19 matrix built in the basis of 
functions constructed from the "left" and the "right" dimers Eq. p^)) 

(StS r ; S lr ||5 cl(dl) || S;s; ; S' lr ) = ^/ (2S lr + l)(2S' lr + l)(-l) 1+s > +s ^ I ^ ^ ^ I (11; S t ||5(1(2))|| 11; S' t ) 5 SrS ,. 

(A2) 

The RME for spins on the c5 (d5) sites amount to 

(5*5V; S lr ||5 c5(d5) || SlfySlr) = ^(2S lr + l){2S[ r + i)(_i)i+s«+*+*«r f ^ J | Sr \\ S {1{2))\\ 11; S' r ) 8 Sl3 >. 

(A3) 

The reduced matrix elements of spin operators on sites c2(c?2), o4(ti4) calculated on the eigenfuctions of the upper 
and down parts form the 73789 x 73789 matrices. 



(^upSupidownSdown; *~' ud 1 1 *"' c2 ( c4 ) 1 1 ^up^up^down^down) S u d) — \J (25 u d + l)(25^ d + 1)( 1) + S "P + S 'down+' S ud 
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Sud 1 ^ud 
^down '-' U P '-'down 



idown'S'down 1 1 <Sc2(c4) 1 1 'down^down ) ^ up ii p ^S„ p S^ p • 



(A4) 



(iup'S'upidown'5'down; Sud ||i5d2(d4)|| ^up'^'up^dowii'^'down! ^ud) _ \j (2S'ud + l)(2S'[ ld + 1) ( — l) 1+Sup+ ' Sdown + ' S,ud 



"Sud 1 ^ud 
'-'up ^down ^up 



^upS'up 11^2(^4)11 ^up'S'up) ^dow»»^ OWI /s down S^ 



(A5) 



The reduced matrix elements of spin operators on sites c2(e?2), c4(g?4) calculated on the eigenfuctions of the envi- 
ronment. The dimension of this matrices determines by dimension of truncated basis of environment 



/'S'cnv H'S'fcll ^env^env) — ^ y A 



X { (^up^upMown^'down) ^ud? (^/^V) "^irj "^env ||^fc|| (Vp^up^down'^dowii) ^ud? (^7*-V) ^Iri ^env) 5 



(A6) 



where /c = c2(d2), c4(d4) and 



own) 'S'udj (SlS r ) Sir] S C nv 1 1 ■S'fc 1 1 ('up^up'down^down) '-'ud I ('-'i'-V) ^lr> ^env) — \Z(2"5cnv ~l~ l)(2S'env -0 



<? 1 f 

^env 1 ^env 

'-'nH ^Jr '-'ud 



(*up5'upidown*5'down; <Sud \\Sk\\ Cp^up^down'Sdown! ^ud) ^SiS'^Sr-S^SirS'^- (A7) 



APPENDIX B 



Quantum statistics is based on two principles, the first is the exchange statistics, when a permutation of two 
identical particles causes an appearance of a phase factor in the total wave function, and the second is the exclusion 
statistics, which reflects an ability to accommodate p particles in the same single- particle quantum state. Whereas 
the first concept depends on the space dimensionality of the system, the second one does not.— 

The exclusion statistics algebra obeying the generalized Pauli exclusion principle can be formulated in terms of the 
bond g operators that has been used in the main text. Another variant of the exclusion statistics can be realized, for 
example, via Green's parafermion statisticsi 21 i 22 According to common formalism based on Burnside's theorem of the 
group theory (see RefJ^ for details) the both algebraic approaches are related with each other. 

Indeed, let us introduce two modes (a = 1, 2) for the each i-th bond 



{df,d«} = {(df)\(d^}=0, 
{d?,(^) f }=% 



with the condition d"\ vacuum) = 0. For a ^ f3 the modes satisfy non-standard relations 



df, (rf; 



= 0. 



Parafermion creation and annihilation operators are determined as 



-;t 



d)+d]. 



They satisfy the commutation relations 



= -25 u d 



d\,d.j ,di 
[[di,dj] ,di] = 0. 



(Bl) 



(B2) 



(B3) 



(B4) 
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The parafermion number operator rij = (dj)^ dj + {df} 1 dj can be written as 

j 1 



2 , 



and obeys the commutation rule 



= %4- 



From the property (»ij) 2 = it follows that nlj varies from to 2. Moreover, 



l\t /j2\t 



2(dJ) T (* 



(B5) 



(B6) 



(B7) 



that means (d^j = 0. Therefore, the parafermion representation provides the extended Pauli exclusion principle. 

To establish a connection between the bond (/-algebra and the parafermion statistics, we note that the local Hilbert 
space related with a bond has the dimension D = 3. Therefore, one can map the (/-particles onto the algebra of 5-1 
operators 



S+ = V2gJ 



V2 



1 nf 



s; 



V2 



1 

71 



1 n? 



(B8) 



Sf =nf-l. 

These spin operators are connected with two-flavour hard-core bosons via the generalization of the Jordan- Wigner 
transformation 23,24 



St = s/2 (bl + b l2 ) , 

Sr = V2(6a + 4). 
S? = blba - b\ 2 b i2 



(B9) 



with the imposed constraint b^bn + b\jbi2 = 1, and the spin state S z = is taken as a vacuum. The commutation 
relations for the hard bosons are 



[bia, b^] 



0. 



bia , 



iff 



= S l3 6 a p (l-n b ia ) 



a ia> u jf3 



= 5i 3 5 a pb\ a , 



(BIO) 



where n\ a = b\ a bi a (a — 1,2) is the number operator for the hard bosons. 

A transition from the hard-core bosons to the parafermions is related with a transmutation of statistics. In 
two-dimensional case the change of statistics is based on a generalization of the conventional Jordan- Wigner 
transformation! 25 ' 26 In the following, for simplicity, we illustrate the connection on an example of dimerized one- 
dimensional S-l chain. 

The parafermion modes are converted into the canonical two- flavour canonical fermions c !Q (a = 1,2) determined 
on the z-th bond of the chain through the partial non-local transmutators 



(4) = 4i ex p 



3<i 



i d i) = c a ex P 




where 



{Ci a ,Cjff} — {ct ,ctg| — 0, jciQDct^l — Sij5 a f3, 



(Bll) 



(B12) 
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and rii a = c ia Ci a is the number operator for the fermions. 

A map between the hard-core bosons and the two-flavour fermions is established by the total non-local transmutators 



The relations (|B8IB9IB11IB13|) provide a map between the parafermions and the g-particles. 
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